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Computational experiments are used to show that grain boundary mobility is independent of 
driving force in a two-dimensional, square-lattice Ising model with Metropolis kinetics. This is 
established over the entire Monte Carlo temperature range. A calibration methodology is then 
introduced which endows the Monte Carlo algorithm with time and length scales and allows the 
Monte Carlo parameters to be expressed directly in terms of experimentally measurable parameters. 
A comparison of results obtained for a variety of driving forces and temperature conditions indicates 
that such Calibrated Monte Carlo models are able to capture the grain boundary kinetics predicted 
by sharp-interface theory. 
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INTRODUCTION 

Meso-scale simulation of polycrystalline grain coarsen- 
ing can be efficiently carried out using Monte Carlo (MC) 
algorithms. Each site represents many millions of atoms 
of like orientation, and this orientation is tracked via a 
discrete, non-conserved order parameter [l|, 01- As op- 
posed to the atomistic kinetic Monte Carlo (KMC) meth- 
ods [3| , meso-scale MC models have no intrinsic time and 
length scales. However, they can be endowed with such 
physicality by developing a parametric correspondence 
with sharp-interface (SI) properties which can be exper- 
imentally measured. The result is a Calibrated Monte 
Carlo (CMC) algorithm [J. 

The calibration takes advantage of the fact that MC 
and SI paradigms both exhibit a linear relation between 
grain boundary velocity, v, and total driving force, P 
Within the deterministic SI setting, this is written 
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with M the grain boundary mobility. The kinetic relation 
can be nondimensionalized using a characteristic length, 
Zq, time, to J and energy, E^. The nondimcnsional kinetic 
quantities, 
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Here and subsequently, parameters without a bar arc un- 
derstood to be nondimensional. In general, driving forces 
are generated by either capillary or bulk effects. Curva- 
ture forces are described by 0, @| 
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where 7 is the surface energy, 7 is the second derivative 
of the surface energy with respect to inclination angle, 
and K, is the grain boundary curvature; (7 + 7 ) is the 
grain boundary stiffness Q, which is also written as 7*. 
The addition of a driving force due to a bulk energy dif- 
ference per unit volume, b, across the boundary gives the 
general kinetic relation 
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This expression is valid provided that grain boundar 
mobility is independent of the types of driving forces 
0| . The experimentally measurable quantities which must 
be supplied to this continuum, sharp-interface model are 
M, b, and 7. These quantities can be used to calibrate 
MC models. 

Now consider a MC approach to study the same prob- 
lem. A two-state field distinguishes grain orientation on 
a fixed lattice with microstructural evolution described 
as a Markov chain [lo| . A nondimensional interfacial 
free energy can be numerically calculated or analytically 
derived. This free energy along with bulk energy differ- 
ences can be used to construct MC driving forces. There 
are many ways in which a spin lattice can be evolved, 
but the Metropolis algorithm is thought to properly 
capture the essence of thermally activated events. A lat- 
tice site is chosen at random and is temporarily given a 
new orientation. The probability of accepting such a trial 
fluctuation of the system is made proportional to an Ar- 
rhenius relation involving the energy change and a non- 
physical Monte Carlo temperature, Tmc- An MC time 
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step, Tmc, can then be defined as a set of N such trials 
with N the number of lattice sites. Computational exper- 
iments can thus be contrived to measure the proportion- 
ality between accretive MC speed and an isolated type of 
driving forces. e. to measure MC mobilities. Such mo- 
bilities will, of course, be a function of MC temperature 
and could, in principal, depend on the driving force as 
well. The MC mobilities may be influenced by the under- 
lying MC lattice and should not be assumed to naturally 
capture the inclination dependence of the physical sys- 
tems they mimic. If the proper inclination dependence is 
not captured, though, a MC model will not produce mor- 
phological evolutions which correspond to experimental 
measurements. 

A significant number of studies conclude that, in fact, 
morphological evolution is correctly captured by MC sim- 
ulations [l|, 0, ■ This suggests that the MC mobilities 
do capture an inclination dependence which matches re- 
ality. The implication is that if a calibration can be made 
for a prescribed boundary geometry then it will apply to 
arbitrary geometries as well. Consideration must there- 
fore be given to the functional form of the mobilities since 
the thermodynamic driving forces can be either derived 
or numerically estimated. 

Recent investigations have made important inroads to- 
wards understanding the mathematical structure of MC 
mobilities and have set the stage for more comprehen- 
sive numerical work. Lobkovsky et al. [1, U, 12| studied 
grain boundary mobility with a two-dimensional Ising 
model on a triangular lattice. The grain boundary 
stiffness, which comprises the capillary driving force, was 
derived under the assumption that the MC temperature 
is sufficiently low that a kink model of the interface is 
valid. To understand the methodology, consider a flat 
boundary with length, L, and inclination angle, (p. The 
number of geometrically necessary kinks is 
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where a\/3/2 is the triangular kink height. Thermally 
generated kinks can be neglected at low temperatures, 
and the total number of lattice sites between the geo- 
metrically necessary kinks is then 
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The resulting configurational entropy of the flat bound- 
ary is 



and the grain boundary stiffness, a*, is therefore 
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This stiffness does not endow the MC model with a time 
scale though. That is accomplished in two pieces. First a 
grain boundary mobility is approximated using the kink 
drift velocity [8| to be 
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Here r is the nondimensional time associated with one 
MC step. The remaining piece of the calibration is to 
prescribe a physical time to each MC step. In practice, 
this time calibration could be accomplished by comparing 
the rate of MC boundary evolution with that measured 
experimentally. From this point, a relationship can be 
established between the dimensional sharp-interface pa- 
rameters and the MC parameters. The resulting CMC 
code would allow experimental parameters to be used di- 
rectly in the MC paradigm and would endow the Markov 
chain with a physical time scale. 

Lobkovsky et al. showed that both the mobility and 
stiffness are anisotropic but that the reduced mobility is 
isotropic This suggests that there may only be a 

single MC mobility to contend with and that its depen- 
dence on inclination angle can be derived from that of 
the grain boundary stiffness. This is consistent with an 
earlier computational study on a 2-D square lattice [l[ 
where the reduced mobility was found to be isotropic as 
well. In a subsequent numerical study at a single (low) 
temperature {Tmc/ J = 0.2), Lobkovsky et al. found the 
mobility to be independent of driving force Q. 

In the current work, we build on these studies within 
the setting of 2-D square lattices. Consistent with 
Lobkovsky et al. Q , we find that the reduced mobility 
is isotropic. It is numerically measured under the action 
of capillary forces, and the mobility is then backed out 
using an inclination-dependent analytical solution for the 
grain boundary stiffness. This mobility is then applied to 
study planar boundaries with a range of inclination an- 
gles that are driven only by bulk energy differences. We 
find that, over a broad range of MC temperatures, the 
mobility is independent of driving force. This extends the 
conclusions of Lobkovsky and co-workers to the square 
lattice system and to the full spectrum of MC temper- 
atures. The results lend additional weight to possible 
universality of two key properties: that reduced mobili- 
ties are isotropic; and, there is only a single mobility that 
is independent of driving force. The assumption is then 
made that both conclusions are valid for experimental 
boundaries in order to focus on calibration methodology. 
Obtaining the characteristic time, the only non-trivial ef- 
fort, is accomplished using a computationally measured 
MC mobility for a single inclination. The work extends 
our previous results by removing the restriction to 
high MC temperatures where isotropy can be assumed. 
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THE SQUARE LATTICE ISING MODEL 

The two-dimensional, square lattice Ising system is 
characterized by an interfacial energy that varies strongly 
with inclination angle-i.e. is anisotropic [14[. For an in- 
terface with zero inclination angle, Onsager derived the 
interfacial energy as a function of temperature [l5| : 



a{a)A 
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(11) 



where a = Tmcl J is the effective temperature, A is the 
side length of one lattice cell, and T^c and J are the 
fundamental temperature and interaction energy, respec- 
tively. This analytical expression for surface energy was 
subsequently extended by Abraham to allow for an arbi- 
trary inclination angle, i^. Til 17. 18. loj: 
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where 



771 = sinh ^ [a((/3, a) cos (/s] 

772 = sinh^"'" [a((/3, a) sin (/j] 
a((p,a) = A[l-(2M)2]i/2 

x{l-h [sin^ 2ip+{2/Afcos^ 2v?]i/2}-i/2 
A = cosh2 2/!r/sinh2i4: 

K = 1/a. (13) 

Figs. [T] and show the dependence of this energy on 
inclination angle and temperature, respectively. Eq. (I12|) 
can then be used to derive grain boundary stiffness of the 
Ising model [20| : 

^ _ (a- + (T )A _ aa{Lp,a) 

J iv3 + m) 

774 = cosV(l + sinVa^)^^^- (14) 

Figs. [3] and m show the grain boundary stiffness, a*, and 
the second derivative of surface energy, a , as functions of 
inclination angle for three values of effective temperature, 
a. At low temperatures, a dominates the stiffness and 
it strongly depends on inclination angle. At high tem- 
peratures, though, the stiffness is almost isotropic since 
(T is negligible. The grain boundary stiffness is then well 
approximated by the surface energy a. 

With an analytical expression for the driving forces 
in hand, the next step is to computationally determine 
the MC mobility function. This is undertaken within a 
purely capillary setting: 



Mj*K = M*K. 



(15) 



A second configuration can then be used to measure the 
speed of a planar boundary driven only by a difference 
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FIG. 1: Temperature dependence of surface energy for three 
inclination angles, ip. 
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FIG. 2: Angle dependence of surface energy at three effective 
temperatures, a. 



in bulk free energy: 



M 6, 



(16) 



The MC mobility can be computed here as well and com- 
pared with the capillary result to determine the degree 
to which the mobility is independent driving force. 
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FIG. 3: Angle dependence of grain boundary stiffness at 
three effective temperatures, a. 
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FIG. 4: Angle dependence of the second derivative of surface 
energy with respect to inchnation angle. This quantity is 
plotted for three values of effective temperature, a. 




FIG. 5: Monte Carlo simulation of inner grain shrinking at 
(a) low temperature without bulk energy, (b) high tempera- 
ture without bulk energy. The grid size is 200 x 200, and 100 
experiments were averaged. The dotted circles are a guide to 
the eye. 



FIG. 6: Numerically derived inner grain shrink rate and bi- 
linear fit. A 200 X 200 grid was used with an initial internal 
area fraction of 0.5. Each point represents the average of 200 
experiments, and error bars show the standard deviation of 
the shrink rate for each temperature. 
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FIG. 7: Illustration of the geometry used to measure the 
mobility of flat grain boundaries driven by an external field. 
The grid size is 400 x 100. 



MONTE CARLO MOBILITY FOR THE 2-D 
SQUARE LATTICE 

Consider a 2D square grain which encloses a circular 
inner grain. In the absence of external fields, SI kinet- 
ics implies that the area of the inner grain A(t), should 
shrink according to 



dA 

— = 2ttAH* 
dt ' 
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(17) 



The reduced mobility therefore corresponds to the shrink 
rate and, once measured, can be used to determine the 
MC mobility. To accomplish this, MC simulations were 
carried out over a wide temperature range (0 < a < 2.0) 
with the first nearest neighbors considered. Consistent 
with results for other lattices Q, the reduced mobil- 
ity was found to be isotropic at all temperatures. This 
is shown in Fig. [5] for two representative temperatures. 
The results of the analysis are plotted in Fig. [51 which 
exhibits a nearly bi-linear dependence of reduced mobil- 
ity on temperature: 
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This is a refinement on a result we previously ob- 
tained The MC mobility therefore is: 



M„(a, cp) 
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with the grain boundary stiffness, a* , is given in Eq. (jl4p . 

Although this mobility was obtained for a particular 
type of driving force, we postulate that it should apply 
to arbitrary driving forces as well. A second bi-crystal 
setting was therefore considered in which a planar grain 
boundary was driven by an external field (Fig. [7]). Two 
parallel, flat interfaces were constructed with an incli- 
nation angle of (p. Periodic boundary conditions (PBCs) 
were enforced at left and right, and screw periodic bound- 
ary conditions (SPBC) were enforced at top and bottom. 
A horizontally oriented external field was applied with a 
magnitude 6,„/A^. The bulk energy density was taken 
to be the product of this field and the local spin value so 
that jump in bulk energy across either grain boundary 
has a magnitude of twice this quantity. A positive value 
of bm therefore causes the two flat interfaces to move to- 
wards each other since the spin value of grain A is -1. 
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FIG. 8: Comparison of grain boundary mobilities measured 
for two types of driving forces. The discrete points were ob- 
tained from computational experiments with planar bound- 
aries while the solid curves were generated using the capillary 
result of Eq. (fT9)) . 



The grain boundary velocity can be easily determined by 
tracking the evolution of grain area, and Eq. (|16p was 
subsequently used to calculate the grain boundary mo- 
bility. Fig. [8] shows the external field driven mobility 
with discrete points, and the capillarity driven mobility 
is shown by solid lines. Furthermore, the convergence of 
measured mobility with the two mechanisms confirms the 
mobility is independent of the type of driving forces. As 
with the surface energy and grain boundary stiffness, the 
mobility becomes isotropic as temperature is increased. 



NUMERICAL UNCERTAINTY IN THE 
MEASURED MOBILITY 



The probabilistic nature of MC simulation calls for 
a careful analysis of numerical convergence associated 
with computationally derived relationships such as the 
reduced mobility. To this end, the capillary shrinking 
of a circular internal grain was considered once again 
to determine the influence of number of experiments, n, 
grid size, L and reduced temperature, a, on the reduced 
mobility. For the purposes of this study, the averaged 
reduced mobility is taken to be a function of these pa- 
rameters, M*{a,n,L). To generate the fit of Eq. 
the reduced mobility was estimated using n = 200 and 
L = 200. The reduced mobility was estimated to be 
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where M* stands for the measured reduced mobility as- 
sociated with the i*'* experiment. The error bars shown 
in Fig. [5] are the standard deviations (S) of the measure- 
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FIG. 9: The affect of the number of experiments on the 
standard deviation of the averaged reduced mobility. 



mcnts at each temperature: 



S{a, n, L) ~ \^ 



1=1 



(M*(a,L)-Ar(a,7i,L))2 



(21) 



This quantifies the statistical character of MC simulation. 
It is then natural to ask if 200 samples are enough to gen- 
erate a mobility function. The issue can be addressed as 
follows. First, choose the number of experiments, n, in 
a data set and compute the averaged reduced mobility 
based on this set of data using Eq. (PD|) . Then, repeat 
this process a large number of sets, m, and calculate the 
standard deviation of the averaged reduced mobility us- 
ing 



- M*''^\a, n, L, m)f 
' m 



(22) 

Here M*{i) is the averaged reduced mobility on the i*'* 
set of data and M*-''^*{a, n, L, m) is the mean of averaged 
reduced mobility of the m data sets: 



M*-'"'*{a,n,L,m) 



(23) 



We carried out a set of 1000 MC simulations at a = 1.0 
with grid size 200 x 200, and treated the set of reduced 
mobilities are viewed as a reservior for statistical analy- 
sis. The number of samples were then tested out up to 
300. The discrete points in Fig. [5] show the standard de- 
viation of sets of averaged reduced mobility with respect 
to the number of samples, and the gray line represents 
the numerical fit provided in the figure. The set-based 
standard deviation demonstrates the convergence of the 
reduced mobility with respect to the number of experi- 
ments performed. 

The impact of grid size on mobility was taken up next. 
A set of 200 experiments were performed with a = 1.0 for 
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FIG. 10: Grid size effect on the standard deviation of inde- 
pendent measured grain boundary reduced mobility. 



a range of grid sizes. The resulting standard deviations 
are plotted as discrete points in Fig. 1 101 with a numerical 
fit shown in gray. This quantifies the degree to which grid 
size can be used to decrease the variation MC simulation. 



PARAMETRIC LINKS BETWEEN MC AND SI 
MODELS AND APPLICATION 

The computationally derived MC mobility can be used 
to establish the parametric links between nondimensional 
SI and MC models. First recall that the the bulk energy 
correspondence is 



25™ 
A2 ' 



(24) 



Within a purely capillary setting, the MC kinetic relation 
can be written as: 



dA AM„ A\ , 



(25) 



Comparison of Eq.s (p5|) and (|17p provides the desired 
link between the nondimensional SI and MC mobilities: 



M 



(26) 



As with the length scales, Tmc is set to one, which simply 
means that one MC step is equal to one nondimensional 
SI time step. Moreover, the grain boundary stiffness is 
expressed in Eq. (jl4p . The links between nondimensional 
length, time, energy, and mobility are thus established. 
To validate these links, MC simulations were performed 
at a high (a = 1.5) and a low (a = 0.5) temperature 
with a combination of bulk and capillary forces, respec- 
tively. Four bulk energies were tested out at each temper- 
ature. The MC parameters are then converted to nondi- 
mensional SI variables using Eq. ([24|) and Eq. ([26|) . 
Therefore, the corresponding SI simulations can be car- 
ried out. At the high temperature, due to the isotropic 



grain boundary properties, grain boundary velocity can 
be written as 

V = M*{a)K + M{a)b, (27) 
and the inner grain area changing rate therefore is 

^ = 2T:M*{a) + 2tt M{a){-f'^b. (28) 

ut TT 

In contrast, at the low temperature, anisotropic grain 
boundary properties must be applied, and the grain 
boundary velocity now is 

V = M {a, ip)^* {a, ip)K + M {a, ip)b. (29) 

The solution to this differential equation was carried out 
numerically over one-quarter of the computational do- 
main so that the vertical position of an interface point can 
be parametrically described by its horizontal position-i.e. 
y = f{x). The coordinates of points on the grain bound- 
ary are described via finite differencing: 



Xi(i + At) = Xi{t) + ViAtcoslipi 
y^{t + At) = y,it) + ViAtsm[ipi] 



(30) 



ipi^ -+ta.n f {x^{t)). 



These are in terms of the local grain boundary normal 
speed, Vi is 

V, = M{T,ip,)-f*{T,ip,)K, + M{T,ip,)b, (31) 
and the discretized inclination angle: 

(32) 

An expression for the discrete curvature is also required: 
= -^"^--(^^ 3. (33) 

The inner grain area can then be calculated as well: 
A{t + At) = ^|y,(< + At)x,+i(t + At)| - 

i 

^|y,(i + At)x,(t + At))| (34) 

i 

The first set of results are shown in Fig. [Til where time 
slices from four separate experiments are given. The solid 
curves are the numerical solution to the SI kinetic equa- 
tion symmetrically reflected from the upper right quad- 
rant into the other three. Agreement between the MC 
and SI predictions is very good. The top two frames are 
for low temperature and show how the bulk energy can 
cause the inner grain to either shrink or grow. These 
same bulk forces also cause the circular grains to grow 
anisotropically. The lower two figures are at a much 
higher temperature and indicate how isotropy is recov- 
ered for this regime. Again, the MC model compares well 
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FIG. 13: High temperature CMC predictions are compared 
with those from the SI model for a range of bulk forces but 
with the same capillary force. The grid size is 200 x 200, 
a = 0.5, and 100 CMC simulations were averaged to obtain 
the results shown. 



FIG. 11: Snapshots of Monte Carlo simulations of inner grain 
shrinking and growing at (a) a — 0.5, and bm ~ 0.02, (b) 
a = 0.5, and 6„ = -0.05, (c) a = 1.5, and b„ = 0.05, (d) 
a = 1.5, and &,„ = -0.05. The grid size is 200 x 200, a = 0.5, 
and 100 CMC simulations were averaged to obtain the results 
shown. 
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FIG. 12: Low temperature CMC predictions are compared 
with those from the SI model for a range of bulk forces but 
with the same capillary force. The grid size is 200 x 200, 
a = 0.5, and 100 CMC simulations were averaged to obtain 
the results shown. 



with the SI model. A more quantitative analysis is pro- 
vided in Figs. ll2l and [T51 where the change in inner grain 
area is plotted versus time for several combinations of 
bulk and capillary forces. Although the parametric links 
between SI and MC can be established at an arbitrary 
MC temperature, high temperature is commonly adopted 
in MC simulation to reduce the anisotropic lattice effect. 
These links can also be applied in polycrystallinc system. 



CALIBRATION OF THE MONTE CARLO 
MODEL 



The established parametric links can be used to cali- 
brate the MC algorithm so that MC parameters can be 
expressed directly in terms of experimental properties. 
The Calibrated MC (CMC) kinetics then have physical 
time and length scales. This is achieved using character- 
istic length, Iq, time, and energy Eq, to nondimcn- 
sionalizc the physical SI kinetics: 



b 
M 



Eo 

MEpto 



(35) 



With the nondimensional correspondence already in 
place, this allow the MC parameters to be expressed di- 
rectly in terms of experimental quantities. We take the 
physical domain size, L to be the characteristic length 
{Iq = L). Then the product of Eq. ^5^2 and Eq. ((SSJla 
gives 



in = 



?gM*(a)A^ 
IF 



(36) 



The characteristic energy, Eq , is derived from the SI grain 
boundary stiffness, a*: 



Eo = 



a* lo 



(37) 



As an example, consider a square aluminum sample with 
a 2 mm edge length at 500 K . A circular internal grain 
sits at its center and occupies half of the domain. The 
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interface is taken to have the following isotropic proper- 
ties (21I: 



M = 10" 



1.0 J/m. 

The reduced mobility, M*, therefore is 

AF = Ma* = 1.0 X lO^^^mVs. 



(38) 
(39) 



(40) 



An external field is applied, which causes a bulk energy 
density difference between the two grains: 



6= 1.68 X 10* J/m^ 



(41) 



(This may be viewed as an idealization of dislocation 
content as well.) The inner grain will shrink due to the 
combined effect of capillarity and bulk energy difference, 
and the process can be studied directly using the CMC 
model. Choose a 200 X 200 grid so that(A = 1/200). Also 
choose an effective MC temperature of a = 1.5. The high 
temperature is chosen to correspond to the assumption of 
isotropy in the experimental data. As mentioned before, 
the characteristic length is 



lo = 2.0 X 10"^ m 
and the characteristic time, to, is 
7\/*(1.5)/2A2 



36 s. 



(42) 



(43) 



which means one MC step is equal to 36 seconds in the 
physical domain. The last value required is the charac- 
teristic energy, Eq: 



Eo = ^^8Ax 10^6 J 



(44) 



The bulk energy density difference, b, is nondimensional- 
ized as 



5 = ^ = 8000, 
^0 



(45) 



which is then transformed to a MC bulk energy term 



6A2 
~2~ 



0.1. 



(46) 



This completes the assignment of MC parameters for this 
Al example. However, all previous analysis were per- 
formed for temperatures sufficiently high that both the 
experimental and MC kinetics are isotropic. In order 
to fully evaluate the degree to which the CMC model 
captures SI behavior, we now take into account the 
anisotropic properties of grain boundaries; that is, the 
same physical grain boundary at a relatively low tem- 
perature is considered. Now both grain boundary energy 
and mobility are anisotropic, which can be written as 



M{(p) and (J*{(p), respectively. Here, the explicit expres- 
sion of grain boundary mobility and energy is absent due 
to the lack of experimental data. But the reduced mo- 
bility. M*, is still isotropic: 



M* 



1.0 X IQ-^^rr^/s. 



(47) 



If a low MC temperature {a = 0.5) CMC model is applied 
to mimic this physical process, the characteristic length 
is still 2 X 10~'^TO. The characteristic time now is 



° " — w — " ' 



(48) 



Therefore, one time step in Fig. [T^is equal to 63 seconds. 



CONCLUSIONS 

A two-dimensional, square lattice, Ising model with 
Metropolis kinetics was used to study grain boundary 
motion in response to bulk and capillary driving forces. 
It was found that boundary mobility is not a function of 
the type of driving force. This is consistent with experi- 
mental measurements as well as analytical considerations 
of triangular lattice systems at low temperatures. The re- 
sults lend credence to the idea that the independence of 
mobility and driving force is a universal characteristic of 
lattice systems. 

The identification of a single mobility function allowed 
us to address a common complaint directed at Monte 
Carlo investigations: they lack a physical time scale. 
Simple geometric settings can be used to calibrate the 
Monte Carlo mobility using experimental data; this en- 
dows the Monte Carlo paradigm with a physical time 
scale. The introduction of more easily obtained length 
and energy scales admits the representation Monte Carlo 
parameters directly in terms of experimentally measur- 
able parameters. This new methodology was extensively 
tested and agreed well with the kinetic results predicted 
by sharp-interface theory. We refer to such physical prob- 
abilistic paradigms as Calibrated Monte Carlo models. 

In the present work, the focus was on a single grain 
boundary within a bi-crystal. However, the resulting 
Calibrated Monte Carlo model can be immediately ap- 
plied to polycrystalline material kinetics. The mobility 
function, and hence the time scale, is the same and only 
misorientation dependent grain boundary energies is re- 
quired. The mobility investigation and calibration pro- 
cedure was restricted to a two-dimensional setting where 
an analytical expression for the grain boundary free en- 
ergy made the analysis particularly straightforward. This 
sets the stage for a consideration of three-dimensional 
systems, though, where no analytical solution for grain 
boundary energy is available and where the finite rough- 
ening temperature make it especially complicated to work 
with grain boundary stiffness 2^, 2^. 
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